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ABSTRACT 

We perform cosmological simulations of galaxies forming at z = 3 using the hydrodynamics grid 
code, Enzo. By selecting the largest galaxies in the volume to correspond to Lyman-break galaxies, 
we construct observational spectra of the HI flux distribution around these objects, as well as column 
densities of CIV and OVI throughout a refined region. We successfully reproduce the most recent 
observations of the mean HI flux in the close vicinity of Lyman-break galaxies but see no evidence 
for the proximity effect in earlier observations. While our galaxies do return metals to the IGM, their 
quantity and volume appears to be somewhat less than observed. We conclude that either wc do not 
adequately resolve galactic winds, or that at least some of the intergalactic metal enrichment is by 
early epoch objects whose mass is smaller than our minimum resolved halo mass. 
Subject headings: intergalactic medium, galaxies: evolution, galaxies: formation, quasars: absorption 
lines, cosmology: theory, cosmology: observations 
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1. INTRODUCTION 

The diverse structure and myriad of elements in the 
intergalactic medium (IGM) are witness to the strong 
role of stellar feedback in galaxies. What is not clear 
is which class of objects is primarily responsible for the 
enrichment. One set of likely candidates are the early 
forming Lyman-break galaxies (LBGs). These galaxies 
are found to be highly clustered, with high star forma- 
tion rates a nd possess strong w i nds w ith velocities up to 
775 km s-i llPettini et alJl20f)ll l20fil . 

Whether LBGs are the cause of the intergalactic metals 
is hotly contested. The issue hangs on whether the met- 
als were produced by massive galaxies such as the LBGs, 
or whether they originated in smaller objects at even 
earlier times, such as dwarf galaxies or population HI 
stars. There i s evidenc e to support both sides. Results 
from [Adclbcrg er et all (2003, 2005) (hereafter AOS and 
AOS) find a cross-correlation between CIV systems and 
galaxies that is very similar to the autocorrelation func- 
tion between LBGs. Additionally, gas that lies within 
40kpc of LBGs contains strong CIV absorption lines 
and there is a similar association with the OVI systems. 
This is evidence that the LBGs might be the source of 
the metal systems. However, studies of the CIV c olumn 
density between th e redshifts of 1.5 — 5 (Songailal l2001t 
iSchave et alJ l2(ffl^ seem to show relatively little vari- 
ation in the amount of carbon in the IGM over time, 
suggesting that the IGM metals were already in place at 
the highest observable redshifts, and must therefore stem 
from the second option of small, early epoch objects. 

One query surrounding LBGs is whether they are able 
to thrust the metals they produce out of their gravita- 
tional pull and into the IGM. Observational studies of 
the LBGs have suggested that outflows from the galax- 
ies are strong enough to bodily displace the gas and 
A03 observed a possible void in the absorbing gas at 
distances out to lh~^Mpc from the galaxy. This prox- 
imity effect results in an increase in the transmitted flux 
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near the galaxy, where one would intuitively expect to 
see a decrease as the galaxy is approached. The most 
likely mechanism proposed for this was strong kinetic 
winds (Croft ct al. 2002), which s eems plausibl e base d 
on the observations of winds by iPettini et al.l \2Q()'3j. 
However, recent simulation work has shown these winds 
have little effect on the absorption at th ese distances 
()Kollmeier et al.ll2005tlBruscoh et al.l2003|) and later ob- 
servational work in this area (A05) has failed to confirm 
the increase in flux. Whether LBGs are therefore capa- 
ble of dispersing metals into the IGM remains an open 
question. 

In this letter we perform cosmological simulations of 
galaxies forming at z = 3. By constructing spectra that 
would be observed along multiple lines of sight from a 
distant quasi-stellar object, we measure the flux close to 
the galaxies and the column densities of HI, CIV and 
OVI in the IGM. By comparing these results with ob- 
servations, we search for signs of a proximity effect sur- 
rounding LBGs and aim to determine the most likely 
cause of metal enrichment in the IGM. 

2. COMPUTATIONAL METHODS 

Our simulations were performed usin g the hydro- 
dynamics adaptive mesh code, Enzo (|Brvanl Il999t 
lO'Shea et ani2004() . We used a box with a comov- 
ing length of 20h~^Mpc and a ACDM cosmology with 
(r!A, riDM, ^^b, /», <tb) = (0.7, 0.26, 0.04, 0.67, 0.9) 

The simulation was evolved through to z = 3 and the 
location of the halos was found us ing a halo finder devel- 
oped by lEisenstein fc Hu3 l)1998|) . Initially, a low reso- 
lution run of the whole box was performed and the halo 
positions located. The simulation was then repeated with 
a section of approximately 5h~^Mpc resolved to a max- 
imum resolution of 1.8 kpc with a dark matter particle 
mass of 5 X 10^ M0. This translates into a minimum 
resolved halo mass of approximately ^ 5 x 10^ Mq. 

The star particles formed in the simulation return ther- 
mal energy and metals back to the baryons. The energy 
is put back into the thermal energy of the gas over a 
dynamical time and is equal to 10~^ of the rest mass 
energy of generated stars, equivalent to one supernovae 
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Fig. 1. — Projections of the baryon density (left) and tlie frac- 
tional metallicity, Pmetal/Ptotal- The white circles on the baryon 
density image show the location of the six Lyman-break galaxies. 
Note that the apparent close proximity of the central two galax- 
ies is misleading, due to the image being projected down the y- 
axis. Image size is 5.31 X 4.69h~^Mpc with an integrated depth of 
4.38h-lMpc. Scales are logarithmic with baryon units MqMpc . 



of 10^^ ergs per 55 Mq formed. It is known that some 
of this energy i s radiated away, a lthough previous work 
with this code ijTassis et al.ll200^ has shown that it can 
eject gas from halos. Metals are returned to the gas over 
the same time period with a yield of 0.02. Radiative 
cooling is computed including metal line cooling based 
on the local metal content. To convert the metallic- 
ity into the ioni zation fractions for CIV and OVI, we 
used CLOUDY IIFerland et aUlTOQl. a ssumi ng ioniza- 
tion spectrums from 'HaarfO£A/Iadau' (199()) and, for 
the softer UV background, Haardt & Madau ( 2001) 

Ideally, the LBGs in our simulation box would be 
chosen based on a calculation of their spectrum in the 
U-band, as with observations. However, for simplic- 
ity we selected the halos whose mass was the same or- 
der as the observed halo masses of LBGs, ioii-5='=0-3 M© 
(A05). With this lower limit, we get a space density 
over the simulation box of 6.75 x 10~^ h^Mpc""^ in good 
agreement with the obse rved value of 8 x 10~^ h^Mpc"'^ 
(|Adelberger et all I1998|) . Within our refined region 
(which is actually overdense with respect to the rest of 
the box), we get 6 halos that are above the cutoff mass. 
These are circled in white in the left-hand panel of Figure 

m 

3. RESULTS 

Figure shows the baryon density (left panel) and 
metallicity of the refined region containing 6 halos we 
classify as LBGs. We can see that the LBGs are situ- 
ated at the intersection of multiple gas filaments, and 
so mark out the densest areas of the simulation box in 
agreement with the obs erved highly biased distribution 
l)Adelberger et al1ll998() . They are also a source of met- 
als, with each region round the LBG being surrounded 
by a high metallicity bubble, as shown in the right-hand 
panel of Figure d 

3.1. Flux properties 

To analyse the distribution of the HI near to the LBGs 
in our simulation, we computed the optical depth along 
2000 randomly chosen lines of sight. By calculating the 
observed distance to each galaxy along the lines of sight. 
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Fig. 2. — Mean transmitted flux as a function of distance from 
the Lyman-break galaxies. The symbols mark out the observa- 
tional results from AOS (open circles) and A05 (filled circles). Our 
simulation results are shown by the solid line and agree well with 
the most recent data set. All plots are normalised to agree with 
the 2005 (A05) observations at = 3.25 Mpc. 



the average flux, (F) = e""^, as a function of distance 
from the galaxy, q, was measured. The flux was divided 
and averaged into just under 50 bins and these results 
are presented in Figure [3 Since the exact amplitude 
of the background ionizing radiation is uncertain, the 
results are scaled to match Adelberger's 2005 results at 
q = 3.25h-iMpc. 

From Figure [3 we see a clear difference in the trends 
of the observational results at low q. The earlier obser- 
vations of AOS (open circles) show a definite increase in 
the transmitted flux, which is missing in the later results. 
In their paper, A05 note that the most likely reasons for 
this are the increase in sample size between AOS and A05, 
or an evolution in the LBGs between the epochs of the 
two data sets, taken at z ~ S and z ~ 2. Our simulation 
results (black line) show no proximity effect like that in 
AOS, but agree closely with A05, indicating that sample 
size was the cause of the discrepancy, but not evolution 
since our data is at the earlier epoch of z ~ S. 

The distribution in absorbing gas can be further stud- 
ied by looking at how it varies within lh~^Mpc of the 
LBGs. Figure El shows the distribution of the flux decre- 
ment, 1 — (F), within this volume for the simulation re- 
sults from Enzo, the recent observations from A05 and 
the simulation results performed by KoUmeier et al., 
(taken from A05). Of the 2000 lines-of-sight that we 
passed through the simulation reflned region, 884 passed 
within 1 h~^Mpc of an LBG. The average flux decrement 
within this area was calculated for each of these lines and 
scaled as for Figure [5] A normalised scale, with the to- 
tal number of measurements the same as Adelberger's 
observations, is marked on the left-hand axis while the 
actual numbers are shown on the right. 

The observations indicate a bimodal distribution, with 
gas either very dense near to the galaxy, or largely ab- 
sent. This, A05 argues, is consistent with the idea that 
anisotropic winds are present in these systems which are 
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Fig. 3. — Comparison between the simulated and observed 
Lyman-a absorption within 1 Mpc of LBGs. The abscissa shows 
the flux decrement 1— < FiMpc >i where < FiMpc > is the average 
flux within 1 Mpc of the galaxy. The left-hand ordinate shows the 
relative (in the case of simulated data) number of measurements 
at a given decrement interval. The right-hand ordinate on the top 
panel shows the actual number of measurement for the Enzo sim- 
ulation. The observational and SPH simulation results were taken 
from A05. 



clearing a path in the gas which roughly half the lines- 
of-sight are passing through. If true, this could explain 
the discrepancy between the results, with the earlier 2003 
observations picking out voids on one side of the galaxy 
that the larger sample set in 2005 averages out. The 
Enzo simulation shows some indication of a bimodal dis- 
tribution, but with very few transparent lines-of-sight. 
In addition, there are a larger number of lines-of-sight 
with very high decrements. KoUmeier's results are sim- 
ilar, also seeing a lack of voids in the gas, but with no 
indication of a bimodal distribution at all. The simula- 
tions are therefore not seeing as strong a population of 
galaxies with the anisotropic voids. 

3.2. Intergalactic metals 

To measure the quantity of metals being ejected into 
the IGM, we calculate the column density of HI, CIV 
and OVI along 500 lines-of-sight. These are plotted in 
Figure 01 along with observational data. We note that 
we do not attempt to model observational effects such as 
confusion between HI and OVI lines, which will lead to 
systematic uncertainties at the low column-density end. 

The top panel of Figure^ shows the column density of 
HI in our simulati on (black line) a longside the observa- 
tional results from iPetitiean et" a,nd IHu et all 
(|p,.995). Our results seem to be systematically slightly 
higher at medium to high column densities and curve 
away at lower values. The former is due to our choice of 
a high density area of the simulation box, where a num- 
ber of LEG candidates were placed. The fall off at low 
densities is most likely due to a limitation in resolution of 
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Fig. 4. — Column density distributions for HI (top panel), CIV 
and OVI (middle panel) and for CIV and OVI calculated with a 
softer UV background along with observational data (as listed). 



the simulation ijBrvan et alJII99^ and low column den- 
sity absorbers (which correspond to small filaments in 
under-dense regions) are underestimated. 

In the middle panel of Figure 0] the column densities 
for CIV and OVI arc plotted with the observational data 
for CIV from ScannaDicco ct al. (2 005j) at z = 2.7 and 
z — 1.9. The match between our results and the obser- 
vations at first seems extremely good, pointing towards 
the objects in the region producing the required metal 
enrichment of the IGM. However, as we have seen from 
the hydrogen distribution above, we have picked an over- 
dense region to study and therefore on average, we would 
expect the metal content to be lower. This implies that 
either the metals are more clustered around the galaxies 
than observed, or not enough metals are being ejected 
into the IGM. 

Softening the UV background, as shown in the bot- 
tom panel of Figure 01 improves things slightly, but not 
enough to match the HI over-density. Further soften- 
ing, or a supersolar yield, could resolve the difference, 
but appear to be unlikely (Aguirrc ct al. 2005). There- 
fore, we conclude that either our galaxies do not have 
powerful enough outflows to drive the metals into the 
IGM, or they are not the main source of IGM enrich- 
ment. If they are not the source, it opens the question 
as to what is causing the metal enrichment. As can be 
seen from Figure ^ contributions to our measured metal 
density come not just from the six LEG halos discussed 
in the previous section, but from all resolved galaxies. 
Therefore, if these are not enough, then smaller objects, 
below our minimum resolved halo mass of approximately 
5 X 10^ Mq, must be responsible. This minimum lies in 
the dwarf galaxy mass scale, so possible objects have to 
be dwarf-sized or smaller. 
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If this result is correct and LBGs are not the source 
of IGM enrichment, then the problem is left as to why 
there seems to be a correlation between LBGs and CIV 
systems (A03, AOS ). A so lution to this was suggested by 
iPorciani fc Madaiil l|20n!Tfl . who propose early- forming 
dwarf galaxies, with the same biased distribution as 
LBGs, as the metal source. The LBGs are gravitation- 
ally drawn to the sa me over-dense sites and form in 
bubbles of old metals. IMadau et alJ l|200H) also find that 
metal enrichment is more likely to come from pregalactic 
outflows than later ones fr om LBGs. Addi ti onally , 
simulations performed by iScannaDieco et all l|20fl,^ 
tried multiple schemes of metal injection but found 
only when the metals were placed in bubbles round the 
galaxies were the observations replicated well. They 
point out that outflows from the smaller gravitational 
potentials from dwarves are a more likely candidate for 
the bubbles than the LBGs. 

An alternative explanation for the LBGs weak effect on 
the IGM is that our feedback prescription is not able to 
correctly reproduce outflows of gas and metals. It is clear 
that we cannot resolve the full multi-phase nature of the 
ISM, although we are currently worki ng to understand 
and i mprove our treatment of feedback ijTasker fc Brvanl 
l200f^^ . A visual inspection of our galaxies reveals several 
outflows, with the maximum velocity in the region being 
750 km s~^, encou ragingly close t o the observed outflows 
of 775 km s"^ bv tPe ttini et alJ l|2002t) . However, only 
two of our six halos showed strong outflows and only in 
one direction, indicating our outflows are neither homo- 
geneous nor ubiquitous. Whether this is a problem is 
hard t o determine : work performed by iKollmeier et al.l 
lj2005|) and Thcun s et all (|2002) on the effect of winds 
on absorption indicates that the lack of strong winds is 
probably not the cause of the decreased flux seen in the 
Figure 13 KoUmcicr et al. ( 2005) found that the gas m 
the IGM is unaffected by the presence of winds and has 
only a small i n ipact on th e optical depth close to LBGs. 
iTheuns et all l|2002 f> and 'B ruscoU et all ("200?) likewise 
found that the winds preferentially expanded into the 
voids, leaving the hydrogen filaments intact. It there- 
fore seems likely that we are correctly modeling the HI 
around the LBGs. 

The lack of metal production is more difficult to judge. 
From Figure^ the metals appear closely centered around 



their source, an effect t hat might be genu ine due to 
pressure from the IGM ijFerrara et alJl20?)5^ or due to 
weak feedback. By using an extremely soft UV back- 
ground and prevent ing heated gas from cooling for 10^ yr, 
iTheuns et al.l |2002) is able to reproduce the observed 
CIV density well, implying that galactic outflows at z = 3 
can enrich the IGM if feedback is mo r e effic ient than pre- 
sented here. However. I Aguirre et alJ lj2005|) find that the 
metals are too inhomogenously distributed compared to 
observations, a symptom that would not change in the 
presence of strong winds. They therefore argue that a 
second, high-z, enrichment mechanism is still needed. 

4. CONCLUSIONS 

We performed detailed simulations of galaxies formed 
at z = 3 using an adaptive-mesh code which included 
recipes for star formation, metal production and SN feed- 
back. From this data, we looked at the transmitted HI 
flux in simulated QSO spectra near large galaxies and 
the production of metals that are returned to the IGM. 

Our calculated HI flux in the neighborhood of the 
galaxies agrees well with the most recent observational 
results (A05) in which the mean transmitted flux de- 
creases near Lyman-break galaxies. While the mean 
properties agree well, there is an indication that the dis- 
tribution of transmitted fluxes does not show the same 
strong bimodality seen in the observations. 

We also generate simulated CIV and OVI absorption 
lines and find that, in our simulation box, the column 
density distributions for these tracers are somewhat lower 
than observed. A visual inspection of the simulations 
shows the metals to be quite concentrated around galax- 
ies (both large and small). We conclude that while we 
do find a significant amount of the IGM polluted with 
metals, the amount falls short of observations, implying 
that either our treatment of feedback does not generate 
sufficiently strong winds, or that objects smaller than our 
minimum resolution of about 5 x 10^ M© are responsible 
for the remaining IGM enrichment. 
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